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Strong lensing statistics in large, z ?C 0.2 surveys: bias in 
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1 INTRODUCTION 

Individual lens modeling has already provided many in- 
teresting constrai nts on the mass profile of clusters 
and galaxies (e.g. |Schecht£r_et_alJ LL997I; |Cj>hjorLalJ 2001; 
[ Rusin fc Mai l200lt [Wvithe et all l2002t iKneib et alJ booa 
llreu fc Koopmansl bo03l : iBroadhurst et alJ I2005T) . In the 
case of galaxy lenses, the measured time-delays between 
the multiple images has been used to determine b ounds 
on the value of the cosmolo gical constant (e.g. iRhed 
Il99ll: iKoopmans fc Trel 12003ft and the Hubble constant 
jKundic et alJ Il997bt ISchechter et alJ Il997l : iFerreras etall 
120051) . The cosmological density parameters, Q.a and Q m 
have also been estimated using a sample of gravitational 



ABSTRACT 

Large current and future surveys, like the Two degree Field survey (2dF), the Sloan 
Digital Sky Survey (SDSS) and the proposed Kilo-Degree Survey (KIDS) are likely to 
provide us with many new strong gravitational lenses. Taking cosmological parameters 
as known, we calculate the expected lensing statistics of the galaxy population in large, 
low-redshift surveys. Galaxies are modeled using realistic, multiple components: a dark 
matter halo, a bulge component and disc. We use semi-analytic models of galaxies 
coupled with dark matter haloes in the Millennium Run to model the lens galaxy 
population. Replicating the selection criteria of the 2dF, we create a mock galaxy 
catalogue. We predict that a fraction of 1.4±0.18 x 10 -3 of radio sources will be lensed 
by galaxies within a survey like the 2dF below z < 0.2. We find that proper inclusion 
of the baryonic component is crucial for calculating lensing statistics - pure dark 
matter haloes produce lensing cross sections several orders of magnitude lower. With 
a simulated sample of lensed radio sources, the predicted lensing galaxy population 
consists mainly of ellipticals (~ 80%) with an average lens velocity dispersion of 164 ± 
3kms , producing typical image separations of ~ 3". The lens galaxy population 
lies on the fundamental plane but its velocity dispersion distribution is shifted to 
higher values compared to all early-type galaxies. We show that magnification bias 
affects lens statistics very strongly and increases the 4:2 image ratio drastically. Taking 
this effect into account, we predict that the ratio of 4:2 image systems is 30 ± 5%, 
consistent with the observed ratio found in the Cosmic Lens All-Sky Survey (CLASS). 
We also find that the population of 4-image lens galaxies differs markedly from the 
population of lens galaxies in 2-image systems. We find that while most lenses tend 
to be ellipticals, galaxies that produce 4 image systems preferentially tend to be lower 
velocity dispersion systems with more pronounced disc components. Our key result is 
the explicit demonstration that the population of lens galaxies differs markedly from 
the galaxy population as a whole: lens galaxies have a higher average luminosity and, 
for a given luminosity, they reside in more massive haloes than the overall sample of 
ellipticals. This bias restricts our ability to infer galaxy evolution parameters from a 
sample of lensing galaxies. 

lens systems (iHelbig et al.| Il999t iMeneehetti et alJ 120051 : 
iDalal et all2005l : IChael2003l) . However, despite great efforts, 
modeling of individual gravitational lens systems has not yet 
provided constraints on any of the cosmological parameters 
comparable to those from other techniques like the stud- 
ies of the Cosmic Microwave Background Radiation. The 
primary reason for this is that there are many degenera- 
cies in the lens models which are still poorly understood. 
Since the details of the mass models have considerable effect 
on the derived cosmological parameters, lens model degen- 
eracies pose a fundamental limitation in exploiting lens ing 
jMoller fc Blainlll998t IWucknitdliool IZhao fc QinlbooaY 

A second possibility to constrain cosmology from lens- 
ing is to use lens statistics. With a large sample of observed 
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lenses, the average mass properties of the galaxy population 
become more important and hence the results are largely 
independent of individual deviations in the lens mass dis- 
tributions. This then also enables constraints on the mass 
properties of the galaxy population as a whole to be ob- 
tained. 

Eve n though cosmology crucially affects lensing statis- 
tics (e.g.lTurner et aill984lFukugita et alll992t iKochanekl 
ll99.4lChaell2003h . several effects complicate the connection. 
The source properties, like source redshift distribution and 
luminosity function, together with magnification bias, all 
affect lensing statistics strongly. There is furthermore, de- 
generacy between lens models, statistical source properties, 
magnification bias and cosmological parameters. Previous 
constraints on cosmological parameters from strong lensing 
are thus not directly comparable in precision to other tech- 
niques like con straints from experi ments such as the WMAP 
space mission JSpereel et alJ 2006). It is perhaps more fruit- 
ful to take the cosmological parameters as "known" quanti- 
ties and use lensing statistics instead to constrain the prop- 
erties of the lensing galaxy population. This is the approach 
that we explore in this work. 

The statistics of strong lensing by galaxies has been 
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studied by several authors in the pa st ||TrrnBr_et_al_ 
i Fukueita et alJll992t FKochanekllT993 iMoller fc Blair 
Chae 2003). Most of these studies calculate a lensing cross 
section analytically using a simple lens model along with 
known mass and luminosity functions. Others are concerned 
with the effect of specific details of the mass models on 
the lens statistics o r properties of the source population 
jHuterer et al.ll2005t) . On larger scales, the arc-statistics of 
galaxy clu sters has been used to constrain cosmological pa- 
rameter s jBart elmann et alJ ll99St iMeneehetti et all l2005t 
iDalal et all 120051) . Perhaps due to the limited sample of 
galaxy lens systems known, the interest in statistical lensing 
on galaxy scales has so far concentrated mainly on predict- 
ing the number of strongly imaged systems using simple lens 
models. However, the number of known lens systems is in- 
creasing steadily and in the near future many large surveys 
are likely to detect a significantly increased number of strong 
lens systems. Such a large number of lens systems could 
potentially constrain masses of galaxies and provide impor- 
tant constraints on the mass evolution of galaxies between 
z ~ 0.1 and z ~ 1. Thus, statistical lensing is potentially a 
very important tool for galaxy evolution studies. 

Detailed statistical lens studies have so far been ham- 
pered by the lack of a clearly defined lens sample. Most 
known lens systems have been discovered serendipitously 
and comprehensive surveys like the Cosmic Lens All- Sky 
Survey (CLASS) are rare. Simple estimates of the lensing 
probability suggest the fraction of sources at z ~ 1 being 
lensed into multiple images is between 10 -4 and 10 -2 . The 
largest current survey, the Sloan Digital Sky Survey (SDSS), 
contains about 500,000 foreground galaxies and about 100 
million background galaxies at higher redshifts, and should 
thus contain at least several thousand lensed galaxies. 

In t his paper, we use semi-analytic models of galaxy for- 
mation l[Kairffjirarm_ej L ^^ ville fc Primackll999l: 
[ Kauffmann et alJ 12004 Bower et alJ I200& ICroton etHl] 
2006) to predict the statistical properties of a large sam- 
ple of strong lenses that would be detected by surveys in 
the radio, optical and infrared wavebands. We concentrate 



here on radio sources that are lensed by low redshift galaxy 
lenses with zi cns < 0.2, as environmental effects are smaller 
for such a population. Two large surveys with similar red- 
shift limits have b een carried out to date, the Sloa n Digital 
Sky Survey CSDSS lAdelman-McCarthv et alJIgOOd) a nd the 
Two Degree Field Survey (2dF IColless et al.lbo'oiTl . Even 
though we use the 2dF as a template survey for this work, 
our results are equally valid for the SDSS. In particular, we 
predict the expected lensing properties of galaxies in these 
surveys. This is done by calculating the lensing properties 
of a mock survey catalogue obtained from the semi-analytic 
models assuming a composite mass model consisting of dark 
matter halo (DM), bulge and disc for lens galaxies. The lens- 
ing cross section, image geometries and magnifications are 
calculated for each galaxy individually using numerical rou- 
tines. From these we then create a sample of lens systems, 
which we then analyse further. 

We describe the use of semi-analytic methods to create 
the lens mass models and outline the method to calculate 
lensing probabilities in detail in §|2| In §|3]we show the pre- 
dicted statistical lensing properties of galaxies in low red- 
shift surveys and the dependence on galaxy properties. In 
§01 we apply specific selection criteria and predict the prop- 
erties of a radio selected sample of lenses with z\ < 0.2. 
Our results are presented in §|SJ including the predictions 
for the lensing rate. In §HJwe summarize our key results and 
discuss their implications. Throughout this paper we use a 
standard ACDM cosmology with Qa = 0.75, Q m = 0.25 and 
H = ISkms' 1 Mpc" 1 . 



2 THE MOCK LENS CATALOGUE 

We use the Mill ennium Run numerical simulation 
JSprineel et aljl 20051 together with semi-analytic modelling 
of galaxy formation to predict the properties of the dark and 
luminous matter in galaxies and thereby derive the expected 
lensing statistics in large surveys at z < 0.2. In the proce- 
dure used for semi-analytic modeling a catalogue of halo 
and substructure is generated at closely spaced redshift in- 
tervals, and for each individual z — DM halo a correspond- 
ing merger tree is i dentified. The prescription described in 
ICroton et alJ i2006T) is used to follow stellar masses, lumi- 
nosities and morphologies of galaxies that assemble in these 
DM haloes. From these, we create a mock low redshift galaxy 
catalogue that is similar to the 2dF North field by applying 
the same selection cuts in magnitude and redshift. We de- 
scribe the creation of this mock 2dF galaxy catalogue below. 
The parent DM haloes are taken from the Millennium Run. 
Selection of simulated haloes is done by defining a backward 
lightcone from z = 0. Even though we use the creation of a 
backward lightcone, we note that for most of the work pre- 
sented here, the creation of a full light cone is not strictly 
necessary, as we treat galaxies as isolated in our lensing anal- 
ysis. In this treatment, we are ignoring the immediate envi- 
ronments of haloes while calculating lensing cross sections 
and image geometries. It is well documented that the pres- 
ence of mass concentrations in th e vicinity of lenses im pacts 
the image separation distribution) Moller & Blain 2001J) and 
we address this interplay in future work. 
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2.1 Making mock observations 

We construct mock observations of our artificial Universe in 
the Millennium simulation by positioning a virtual observer 
in the simulation at zero redshift and finding those galaxies 
which lie on the observer's backward lightcone. The main 
limitation in producing a mock observation of a simulation 
is the finite box size, which in this case is 500 h~ Mpc on 
a side. The periodic nature of the simulation allows us to 
fill the space with any number of boxes we require but we 
still have to deal with periodic replications that would ap- 
pear if we were to look through the simulation volume along 
one of its preferred axes. We can avoid this kaleidoscopic 
effect through slanting the observed cone by a certain angle. 
After having determined the survey geometry, observer posi- 
tion and line-of-sight we fill the four-dimensional Euclidean 
space-time with a grid of simulation boxes. In the three spa- 
tial coordinates we make use of the periodicity of the simula- 
tion whereas the time coordinate is given by the 64 snapshot 
times corresponding to their respective output redshifts. In 
practice, only those cells in the space-time grid are popu- 
lated with galaxies which actually intersect the backward 
lightcone in the observed field-of-view. After coarsely filling 
the volume around the observed lightcone with simulation 
boxes in this way one can simply "chisel off" the protruding 
volume. 

We note that by using comoving coordinates and as- 
suming a flat Universe we have the luxury of cutting the 
lightcone from our simulated volume simply like we would 
in Euclidean geometry. In general this is a non-trivial en- 
deavour that requires taking into account the curvature of 
the Universe as well as its expansion with time. 



2.2 Lensing probabilities for an SIS model 

Before calculating the lensing cross section for realistic 
mass models of lensing galaxies we first discuss the sin- 
gular isothermal sphere (SIS) model that has in the past 



most commonly been used 


for statistical lensing 


tions (e.g. iTurner et al] 


ll984UFukueita et al.lll992l) 



its simplicity the SIS model has been surprisingly success- 
ful in describing the overa ll statistical propert i es of galaxy - 
scale lenses. As shown by iKoopmans fc Treul J2002L 12003 ) . 
the stellar component of lens galaxies has a density profile 
that is much steeper at than the dark matter at small radii 
so that the total mass profile (including dark matter) fol- 
lows an isothermal profile on scales probed by galaxy strong 
lensing more closely, thereby explaining the modelling suc- 
cess of the SIS. For most groups and cluster lenses, how- 
ever, the SIS model is an inadequate description of the total 
mass distribution. On cluster scales the baryonic compo- 
nent becomes less important and the lensing profile is more 
accurately modeled using a N avarro, Frenk & White pro- 
file fNFWlNavarro et al.lll997li well beyond the arc radius 
jKneib et al.ll200d) . For lower redshift lenses, the Einstein 
radius is small compared to the galaxy size and the rela- 
tive DM mass content within the Einstein radius is smaller. 
Thus, for all galaxies (excluding the most massive groups 
and clusters), the baryonic galaxy component of the mass 
profile dominates the lensing cross section at low redshifts. 

The lensing cross sections of a SIS model for elliptical 
galaxies can be calculated directly from the observed Faber- 
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Figure 1. The lensing cross section for SIS lenses with sources at 
2 = 2 for different slopes of the luminosity function and power-law 
indices relating luminosity and velocity dispersion. The solid line 
shows the results for a = —1.21 and j3 = 0.25, corresponding to 
the 2dF results. The dashed, red line is for a = —2 and j3 = 0.25 
and the dotted green and dot-dashed blue lines are for a = —1.21 
and f3 = 0.31 and j3 = 0.5 respectively. The functional form of 
the luminosity function is given in equation 151 



Jackson relation. The lensing cross section of an SIS is given 
by: 



TSIS 



(1) 



the total lensing probability can be calculated using the lu- 
minosity function N(L) of galaxies: 



roc / a(L) N(L) dL, 
Jo 



where 



a(L) 



(2) 



(3) 



is the Faber- Jackson relation . Using the 2dF luminosity 
function iNorberg et alJ i2002T) fit, we have: 
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where <j>, = 1.64 x 10" 2 Mpc" 3 , a = -1.21, L, = 1.82 x 
1O 1O L and r, = (a./kms" 1 ) 4 x 1.055 x 10~ 9 arcsec 2 is the 
characteristic lensing cross section of an L* galaxy for equal 
distances Dls and Dos- The comoving volume of a shell of 
thickness dz is given by dV co /dz x dz. 

We show the t — L relation for various choices of a 
and (3 in Fig.0 For a slope of the Faber-Jackson relation 
of j3 < 0.3 the lensing cross section t(L) is a monotonically 
decreasing function with L. In the case of an SIS lens, the 
largest contribution to the total lensing cross section is ex- 
pected from the galaxy population around L t . Using the 
Faber-Jackson re lation determined fro m early-type galax- 
ies in the SDSS feernardi et alJl2003all . with /3 = 1/4 and 
<t„ = <t{L*) — 166.44 kms" 1 , we obtain a total lensing cross 
section of: 



Ttot = 



t(L) dL = 8.36 x 10" 



(5) 



4 M oiler, Kitzbichler, Natarajan 



for sources at z SOU rces = 2. This is our first and simplest esti- 
mate of the expected lensing rate. At this juncture, we have 
neglected the fact that galaxies contain DM and baryons 
which are, in general, not expected to be distributed isother- 
mally. The measured velocity dispersion may be biased with 
respect to the velocity dispersion inferred for the dark mat- 
ter. In addition, we have also neglected magnification bias 
and the source luminosity and redshift distributions. 



2.3 Towards realistic lens models: dark matter 
and baryons 

More sophisticated mode ls than the SIS have been used 
by e.g. lOeuri et al.l (120051) who predicted the lensing statis- 
tics in the SDSS using a spherical Hernquist model for the 
baryons and an NFW dark matter mass profile. Spherical 
models do predic t the overall lensing c ross section with ad- 
equate accuracy ijFukugita et al.l lT992'). They cannot, how- 
ever, predict any detailed properties, like the average magni- 
fication and the statistics of image geometries. For instance, 
the number of four image systems predicted by spherical 
lens models is zero. We use the ray-tracing code gLens 
( Mollcr & Blai nll99ct) 1 to calculate the lensing properties of 
multi-component lens systems efficiently. In order to calcu- 
late expected lensing properties of galaxies in large surveys, 
we need to relate the observed properties, like luminosity 
and size to a mass model. Since lensing is sensitive to the 
total mass, we need to include the contributions both of the 
baryonic and of the dark matter. Our galaxy model there- 
fore consists of a baryonic component, modeled as a bulge 
of mass Mb plus a disc of mass Md that reside in a DM halo 
of virial mass M v i r . 

In this paper, we consider only the lensing statis- 
tics of isolated lenses. In reality, most lens galax- 
ies are expected to reside in dense environments, in 
groups or clusters <|Kundic et al1ll997al iMoller et alJl2002t 
iKeeton fc Zabludofi1l2004l: lOguri et al.ll2005|) . However, the 
lens galaxies we consider here are all at relatively low red- 
shifts, where the effect of environment is expected to be less 
strong. At redshifts z ^ 0.2 the expected Einstein radius 
is small compared with the distance to nearest neighbours 
(this is true even for lens galaxies in groups and clusters). 
We note that by construction our catalogue also includes 
the total lensing cross section by groups and clusters: the 
haloes of these objects, together with their central galaxy, 
are included in our catalogue as single, isolated galaxies. In 
this sense, our calculations are complete, except that they do 
not take into account the effects of satellite galaxies and ex- 
ternal shear deriving from neigbouring mass concentrations 
and halo substructure within these DM haloes. 



2.3.1 Modeling the dark matter halo 

N-body simulations of structure formation in a cold dark 
matter dominated Universe find that the best-fit density 
profile for dark matter haloes follows a universal functional 



1 The gLens software was developed by Moller to solve the lens 
equation on an adaptive grid on the source plane. This ray-tracing 
code permits the addition of an arbitrary number of mass com- 
ponents for lenses. 



form over a broad range of scales jNavarro et alj [l997') . The 
mass distribution for a spherical NFW halo is given by, 

S c pc(z) 



p(r 



1 + 



(6) 



where r s — R vlv /c is the scale length, p c (z) is the critical 
density for closure at the galaxy redshift and 



8c = 



200/3 c 3 



— • (7) 

log (1 + c) - c/(l + c) y ' 

The concentration parameter c is ultimately a unique prop- 
erty for each halo in simulations, however correlations (with 
some scatter) between the concentration parameter and 
the mass or virial velocity of the halo have been fou nd 
jNavarro et al.ll99llB"ullock et all200ll : lEke et al.l200 ill . A 
closer investigation of the strong lensing by NFW haloes re- 
veals that the lensing cross section depends crucially on the 
concentration parameter. For the same halo mass, a halo 
with c = 20 has a lensing cross section that is ~10 times 
higher than a halo with c = 10, and ~ 100 times higher 
than a halo with c = 5. This implies that the scatter in 
these correlations between halo mass and central concen- 
tration parameter is very important for lensing - most of 
the lensing cross section contribution derives from massive, 
high concentration 'outliers'. Due to this effect, the lensing 
statistics of a galaxy population can only be correctly pre- 
dicted if the full scatter in the distribution function of the 
concentration parameter is taken into account. In this work, 
we calculate the concentration para meter explicitly for e ach 



we calculate tne concentration parameter explicitly tor eac 
halo. For this, we use the equation iNavarro et al.l[l997fl . 



where 

A c = ln(l + c) 



1 



(8) 



(9) 



For each halo in the simulation, the parameters V max and 
V v i T are known and we use the eqns. (8) and (9) to calculate 
the concentration parameter for each halo. This procedure 
ensures that the distribution of concentration parameters 
automatically has the correct scatter. 



2.3.2 Modeling the baryonic components 

The baryonic mass is distributed into a bulge and a disc 
component. The relation between the mass of a DM halo and 
the luminosity of the baryonic component can be studied 
using semi-analytic mo dels of galaxy f ormat ion. The semi- 
analytic prescription of ICroton et all {2006) gives a bulge 
luminosity, a disc luminosity and a disc size. 

We show the mass-luminosity relation for our mock cat- 
alogue in Fig.|2| For comparison, we overplot the relation 



L B] /W 10 L Q = 0.3 



(Mvir/lO 1 ^©) 4 



(0.57 + (Af vir /1O 11 M )O-79) 1 /O.23 ' 



(10) 



derived bv lVale fc Ostrikeil (120041) as the solid line, including 
the scatter (dotted, blue line). Note that there are galaxies 
in the catalogue that consistently fall above and to the left 
of the analytic relation (given in eqn. 10) for high virial 
masses, which correspond to central galaxies in the more 
massive haloes of clusters or groups. 
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Figure 2. The relation between the halo virial mass M v i r and 
galaxy luminosity in the Bj-band. The shading shows the num- 
ber of galaxies i n the mock 2dF ca talogue created from the 
Millennium Run ISpringel et al.ll2005t) combined with the semi- 
analytic model of Tcrotori^t^l!n200ri) . The black line show s 
the relation in the Bj-band obtained bv lVale &: Ostrikerl 12004) . 
The dotted blue lines show the 68, 90 and 99 percent levels of 
the scatter around th e mass-luminosity relation as obtained by 
ICoorav fc Cenl I2005T) . 

2.3.3 The Bulge component 

The size and velocity dispersion of the bulge component 
is not determined in the semi-analytic modelling. How- 
ever, the distribution of bulge sizes is known observationally 
ide Jong et alJl2004) . There are no models that relate the 
bulge size to properties of the DM halo. The fundamental 
plane relation suggests that DM haloes and the bulge prop- 
erties are related: virialised haloes with constant M/L will 
naturally lie on the fundamental plane, but so far no the- 
oretical model has been developed that naturally explains 
both the existence and the slope of the fundamental plane. 
In the absence of such a model, we rely on the observed rela- 
tions between luminosities and sizes of elliptical galaxies and 
apply these to the population of bulges as a whole, irrespec- 
tive of the morphologi cal type of the galaxy. In the SDSS, 
iBernardi et, alJ (|2003ar ) find a correlation between half-light 
radius, and galaxy luminosity of the form: 

log-Rb — 0.52 = — log (L/L«), (11) 
a 

where a = 1.5 ± 0.06 and L, = 1.58 x 10 10 L©. They also 
find the Faber- Jackson relation for the velocity dispersion 

log cr- 2.197 = /? log L/L,, (12) 

where j3 = 0.25 ± 0.06. Since the semi-analytic model fixes 
the total bulge luminosity and DM properties, we do not 
use the latter relation. Instead, we only use the Rb — L 
relation (eqn. 11 above) to determine the bulge-size and 
then solve the Jeans-equation to obtain the line of sight 
velocity dispersion a, assuming a constant mass-to-light ra- 
tio of F = 7 in the B-band, which is the m ean M/L ratio 
expected for bulges and ear ly-type galaxies iGerhard et alJ 
l200ll: ICappellari et aljEo 06 ) . This choice does not affect the 
tightness of the fundamental plane relation strongly; the un- 
certainty in q is the strongest contributor to the scatter in 
the relation. We can then check whether the bulges still sat- 
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Figure 3. The Faber- Jackson relation: the number of galaxies in 
the mock catalogue as a function of L and a are shown, with the 
blue line showing the best-fit observed relation a = a-*(L/L*) ' 25 . 
The inset shows the distribution of velocity dispersions, log in 
our mock 2df catalogue. The smooth line in the inset shows the 
observed velocity disp ersion distribution in the SDSS found by 
IBernardi et alJ <2003rl . 
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Figure 4. The fundamental plane for galaxies in the mock cata- 
logue: the number of galaxies in the mock as a function of and 
a combination of velocity dispersion and surface brightness. The 
blue contours show the same for the sample of elliptical galaxies 
in the SDSS iBernardi et al.ll20n.'i3) . The histogram in the inset 
shows the distribution of log in the mock catalogue. The dots 
show the mock lens sample created in §0 

isfy the observed a — L relation and fall on the fundamental 
plane consistent with observed bulge properties. 

We show the Faber- Jackson relation and the histograms 
of the distributions of the velocity dispersion and bulge size 
in Figs.[3]&[I]for the mock catalogue. These plots show that, 
using only the observed R^ — L relation as input, we recover 
the observed Faber-Jackson relation and the fundamental 
plane. Of particular importance is the fact that we recover 
the observed velocity dispersion distribution measured for 
early-type galaxies by the SDSS. We conclude therefore that 
our mock galaxies are realistic. 

These structural parameters are then used to model 
the bulge component in the lens as a de Vaucouleur pro- 
file (Ide Vaucouleurslll948l) . 
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Figure 5. The distribution of disc sizes in the mock 2dF cat- 
alogue. The histograms show the distribution of logiJd in the 
mock, the thick, red is for all galaxies and thin, black is for spiral 
galaxies only. Note that both histograms have been normalised 
independently. 
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Figure 6. The distribution of number of galaxies as a func- 
tion of Vmax and disc luminosity in the mock catalogue. The 
dashed line shows the ob served Tully-Fisher r e lation, Mr. = 
-4.125 - 7.27 x log V m ax iTullv fc Fisherl Il977t iPierce fc TulTvl 
Il992t f 



E b = E ,b x exp-(r/r s ) 1/4 , (13) 

where r s = 3 681 x R b x 0.551 and E 0jb = FL/8n x 
5040rg vl — e 2 , measured in kpc and M^kpc -2 respectively, 
and where F is the fixed mass-to-light ratio of bulges taken 
to be 7Mq/Lq in the B-band. We also allow for ellipticity, 
e, in the bulge. To assign ellipticities to the bulge compo- 
nent in our mock galaxy catalogue, we draw an ellipticity 
randomly from a Gaussian distribution of mean e = 0.3 and 
standard deviation, cr e = 0.3, co nsistent with the observe d 
ellipticities of early-type galaxies llJorgensen fc Franxl l994'l. 



2.3-4 Disc component 

In the semi-analytic models of galaxy formation every 
galaxy, irrespective of morphology, also contains a disc com- 
ponent. The infall and merger history of haloes determines 
the disc properties: total disc mass, M<j = TLd and disc size, 
Rd- The mass distribution of discs is taken to be well- fit by 
an exponential: 

E d = E exp(-r/i? d ), (14) 

where Eo = Md/itR^ is the central surface mass density of 
the disc. 

We plot the distributions of the structural disc param- 
eters for the mock catalogue in Fig.|S] Identifying the maxi- 
mum velocity dispersion Vmax of the DM halo with the disc 
rotational velocity, v c also gives the Tully-Fisher relation for 
galaxies in our mock catalogue (see Fig. 6). 

The blue li ne in Fig.El shows th e Tu lly-Fisher rela- 
tion o bserved bv lTullv fc Pierce! feOOOft and lPierce fc Tulivl 
(1992). There is a very good match between observed and 
predicted relation for Mbj > —21, corresponding to L < 
2.15 x lO 10 !/©. At higher luminosities, there is an upturn 
in the predicted relation. Given the small observed number 
of objects with measured rotational velocities at these lu- 
minosities this is consistent with observational data at the 
present time. 
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Figure 8. The lensing cross section as a function of magnifi- 
cation /i for DM, DM+bulge and DM+bulge+disc models of a 
typical L = lO 11 !;© galaxy in the mock catalogue. Cumulative 
histograms of the cross sectional area for lensing of point sources 
by magnifications of more than fio, t(/i > fio) are shown. The 
solid, red line is for the pure DM model. The solid, black line 
shows the result of adding a bulge and the solid blue line is for 
the complete DM+bulge+disc model. The dashed black and blue 
lines show the magnification cross sections for model of bulge and 
bulge+disc without dark matter, respectively. Model parameters 
are as for the previous figure. 



2.4 Lensing properties of the composite model 

To illustrate the general lensing properties of the composite 
mass model used, we show the image and source magnifica- 
tion maps of the dark matter+bulge+disc models in Fig.|7| 
for a spiral galaxy with L = 10 L§ . Note the strong asym- 
metry introduced by the disc component. This is expected 
since for inclination angles of 9 > 65 deg, the projected sur- 
face mass den sity along the major a xis is strongly increased 
for thin discs jMoller fc BlairJI 1998ft . 

The magnification curves in the source and image planes 
are shown in Fig.|H| The NFW dark matter profile on its own 
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Figure 7. The image (left) and source (right) magnification maps for dark matter+bulge+disc model of a luminous L = 1.8 X 10 11 Lq 
galaxy. The colour scale shows the magnification of a point source as a function of image and source position respectively. The source 
plane is at z B = 1 and the image plane is at z = 0.15. The lens has a bulge to total mass ratio B/T of 0.2, a bulge ellipticity of e = 0.4 
and a disc inclination angle of 8 = 70 deg. The dark halo mass is M V \ T = 1.3 X 1O 13 M0 inside R v i T = 455 kpc, with a concentration 
parameter of c = 9. 



does not have a significant lensing cross section in this ex- 
ample, as the concentration parameter c of the dark matter 
halo is too low. It is interesting to note, however, that de- 
spite this the magnification cross sections for the composite 
model are influenced significantly by the DM component. 
This can be seen by comparing the lensing cross-sections as 
a function of magnification shown as solid (DM) and dashed 
(no DM) lines in Fig.|S] The cross section for the case with 
no DM is a factor of ~ 2 below that with DM, for almost 
all magnifications above about 3. Therefore, even in cases 
when the DM halo is not concentrated enough to produce 
lensing, it still affects the lensing cross section significantly. 



3 THE EXPECTED LENS PROPERTIES OF 
THE GALAXY POPULATION IN THE 2DF 

3.1 Calculation of lensing properties 

Once each galaxy has been assigned a NFW 
halo+bulge+disc mass model, the calculation of the 
lensing probability for a given source population depends 
only on the geometric parameters, i.e. the redshifts of the 
source, lens and the cosmological parameters. We are not 
interested in constraining cosmological parameters with 
lensing in this work and we fix the underlying cosmology as 
described in §Q 

We calculate the lensing properties of each of the ~ 
70, 000 galaxies in the mock catalogue using ray-tracing. 
Since we consider only one lens plane, this reduces to nu- 
merically solving the lens equation, 

9 = j3 -a, (15) 

where 9 is the imaged position of a source at angular position 
(3 and a is the deflection angle. 

For the calculations in this section, we fix the sources 
at Zsourcc = 1. In order to allow ray-tracing calculations of 
sufficient accuracy in a reasonable amount of time we use 
the ad aptive grid technique, as described in Moll er fc Blainl 
teOOll) . For each galaxy, an approximate lensing region size 
is determined using the extent of the Einstein radius. This 
is done by finding the points in the image plane at which 



the magnification diverges. The physical size of the puta- 
tive lensing region is then set to 20 times the distance of 
the most distant such points to the galaxy, giving a conser- 
vative upper estimate of the region within which multiple 
imaging can occur. This region in the image plane is then 
covered with an initial grid of 40 x 40 rectangles. Each grid 
point is then mapped to the source plane. After mapping, 
magnifications on the image plane are given by the ratio of 
the areas of the original to that of the mapped rectangles. 
Depending on the magnification \i the initial grid is then 
adaptively refined by subdividing each rectangle into n su b 
smaller rectangles, with n su b(/x) = 1 -+- {p, — 5)/195 x 20, 
for jj, > 5, and with no refinement otherwise. This resulting 
adaptive grid is then also mapped onto the source plane. Fi- 
nally, the image magnifications and image multiplicities are 
calculated on the resulting fine grid on the source plane. 

3.2 Total lensing cross sections and image 
separation distributions 

Two statistics have b een used in the past to cons train cos- 
mological parameters Jchac 2003; Ko^hanelj[l993j) and len s 
galaxy population properties dChaell2005l : lKochaneklll996l) : 
total lensing cross sections and image separation distribu- 
tions. The total lensing cross section we obtain from our 
mock catalogue is 100752 arcsec 2 in a total survey area of 
657 deg 2 , with 79042 galaxies below redshift of 0.2 consti- 
tuting our lens galaxy population. This gives a total proba- 
bility of a source being lensed as 1.13 x 10 -5 . This estimate 
neglects magnification bias and we show in §2] that magni- 
fication bias increases the lensing probability dramatically. 

We plot the lensing cross section for our mock catalogue 
as a function of lens galaxy luminosity and halo virial mass 
in Fig.|5| The contours show the lensing cross section for 
background sources at z SOU rce = 1 as a function of mass, 
M v i r and luminosity Lb • Comparison with the underlying 
distribution for all galaxies in the sample shows that the con- 
tribution to strong lensing cross section comes preferentially 
from massive, luminous galaxies. The shift in the luminos- 
ity from a peak luminosity of Lj^ t ak ~ 10 10 Lq for the total 
distribution to a peak of = 2 x 1O 1O L0 for the lensing 
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Figure 9. The lensing cross section as function of luminosity and 
mass. The colour scale shows the distribution of all galaxies as a 
function of L and M v i v . The blue contour shows the lensing cross 
section distribution for the complete DM+bulge+disc models. 
The thick, black contour shows the cross section for DM+bulge 
models and the thin black contours are for DM only models. The 
contour levels are 100, 300, 700 and 1000 arcsec 2 for the first two 
cases and 0.1, 1 and 10 arcsec 2 in the latter case. 
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Figure 10. The expected image separations in the 2dF. As in the 
previous figure, we show the results for the DM only (red), the 
DM+bulge but no disc (black) and the DM+bulge+disc (blue) 
models. The solid lines correspond to the image separation his- 
togram of 2 image systems, the dotted line shows the maximum 
image separation for 4 image systems. In all cases, shown above 
the effect of magnification bias is not included. 

sample is clearly strong. An additional shift to higher virial 
masses by a factor of ~ 2 is also evident. The baryonic com- 
ponent dominates the lensing cross section of galaxies, but, 
for a fixed luminosity, the contribution of high mass objects 
to the lensing cross section is larger. The dark matter com- 
ponent therefore also plays a key role. In Fig. llOl we show the 
expected distribution of image separations. Here, we clearly 
see the relative effects of the various lens components. For 
small image separations (corresponding to low masses) the 
cross section is completely dominated by the baryonic com- 
ponent. Large image separation systems, mainly correspond- 
ing to groups and clusters, have lensing cross sections that 
depend much more on the dark matter content. 
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Figure 11. The average ratio of the lensing cross section for 4: 2 
image-systems as a function of image separation. The black line is 
for a DM+bulge only model for all galaxies in the mock catalogue, 
and the blue line for a DM+bulge+disc model. Magnification bias 
is not included in this calculation. 

It is clear from these plots that the baryonic compo- 
nent is absolutely crucial in determining statistical lensing 
properties on galaxy scales. The importance of the baryonic 
contribution to the lensing cross sections of DM haloes is 
clearly seen: it increases by a factor of 10 — 1000 for systems 
with image separations between ~ 0.1 and ~ 1". The in- 
crease is smaller for systems with large image separations. 
These systems correspond to lensing by massive central clus- 
ter and group galaxies. As the Einstein radius increases for 
more massive systems, the fraction of DM within the Ein- 
stein radius also increases, due to the shallower slope of the 
DM component. Despite this importance of the baryonic 
component, the lensing cross section at a fixed L does de- 
pend on the DM halo mass, M v i r , even for highly luminous 
galaxies (cf. Fig.0. 

3.3 Image geometries: two and four image systems 

Previo us work by e.g. Moller fc Blai r] il99ST) and iBlainl 

( 1998) has shown that baryonic components can significantly 
change the total lensing cross section. The disc component 
plays a major role in this; inclined discs can increase the 
lensing cross section significantly. However, as suggested by 
Moll er et alJ i2003l) this effect is important only for late- 
type systems, whereas early-type galaxies have too small a 
disc component to change the lensing cross section signifi- 
cantly. The ratio of 4 to 2 image systems 2 , however, may 
be affected strongly, as even relatively low mass discs can 
increase the asymmetry of the lensing potential noticeably. 

We show the ratio of 4 to 2 image systems in the cat- 
alogue for the DM+bulge and the DM+bulge+disc models 
in Fig.mi For systems with lO 11 /,© < L < 10 12 L o the frac- 
tion of late-type galaxies is high enough to increase the cross 

2 Note that strong lensing produces an odd number of images, 
one of the images is always strongly de-magnified, therefore 2 
images systems have a total of 3 images with 2 magnified ones 
and one central demagnified one and likewise a 4 image system 
has in fact a total of 5 images. 
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section for 4-image systems by more than a factor of ~ 10. 
Again, this result neglects the effect of magnification bias. 
Magnification bias is expected to favour late-type galaxies 
as 4-image lens systems, given that the expected cross sec- 
tion for magnification of fj, ~ 10 — 20 is a factor of up t o ~ 10 
larger for inclined spiral lenses jMoller fc Brainlll99Sft . How- 
ever, as we show in §|1J this is not the case: elliptical galaxies 
have a smaller cross-sectional area for magnifications in that 
range, but they generally have a larger cross section for mag- 
nifications above 20. This implies that magnification bias in 
fact favours the moderately elliptical gravitational potential 
of early- type galaxies. 

3.4 Selection biases 

The results presented above do not include the effects of 
selection biases. All surveys are by construction limited in 
flux, resolution or some other observational criteria. These 
varied selections can strongly affect the inferred properties 
of lens samples. A flux limited survey for example, will be 
strongly affected by magnification bias for sources with a 
steep luminosity function. Since galaxies and radio sources 
follow quite steep luminosity functions magnification bias is 
important for lensing in the radio as well as in the optical. 
One can estimate the effect of magnification bias for sim- 
ple selection criteria. The cumulative probability distribu- 
tion function for magnification above a value A is generally 
P(n > A) oc A~ a , where a — 2 for SIS lenses. Assuming 
that the probability that a source enters the survey is 1 for 
(i x F > Fo and otherwise, it follows that the fraction 
of sources of luminosity F that are obs erved at flu x F' is 
N( F') = f °° dP/dA ( A)N( F' /A)dA [cf. iMaoz et al.1 (Il993h 
and iKing fc Brownd (11996ft ]. For a lens population with a 
more complicated lens model, an analytic calculation of the 
effect of magnification bias becomes infeasible. Including the 
selection effects due to finite resolution limits requires a nu- 
merical calculation - which we describe in the next section. 



4 A MOCK SAMPLE OF LENS SYSTEMS 

We use the lensing code gLens together with our mock 
galaxy catalogue to create a sample of lens systems, includ- 
ing observational selection effects. We include two key selec- 
tion effects: a flux limit and a resolution limit. Concentrating 
on lensed radio sources, we create a sample that has selec- 
tion criteria in concordance with the CLASS survey JChad 
12003ft : 

(i) The total flux F of each system of images is above Fo, 

(ii) the image separation is > 0.03", 

(iii) and for double image systems the magnification ratio 
of bright to faint image is > 0.1. 

We create the sample using rough estimates of the lensing 
cross section for each galaxy assuming spherical symmetry, 
by solving the lens equation for the approximate Einstein 
radius, 9 eT . We then populate the total estimated cross sec- 
tion area, Attempt with Attempt sources. For each source, 
we pick a random flux from a distribution N oc F~ 21 and 
assign a source redshift from a redshift distribution that is 
a Gaussian with mean z s = 1.27 and width 0.95 as found 



by IWillott et ail fcOOlft in the 6CE and 7CRS radio sam- 
ples. For each source, we find the positions of the images 
using an a daptive lensing t echnique similar to the one de- 
scribed in IWucknita (|2004j). We then apply the selection 
criteria (i)-(iii) to the lens system. If the system satisfies 
all these criteria, it is included in the lens sample, other- 
wise another foreground galaxy and source pair are selected 
randomly until the desired number of lens systems is ob- 
tained. For a number of Attempt 'trial' sources, placed in 
an area Attempt << A survcy where A survoy is the total area 
of the survey, there is a number of A^niensed sources that 
are unlensed but above the flux limit for detection. The 
number of expected lenses in the survey is then given by 

Aliens = /lensA^unlensed X A surve y /A a ttemps , with f\ens being 

the lensing fraction. 



4.1 Properties of the lens sample 

We use the procedure described above, to create a sam- 
ple of Anions = 400 lens systems from our mock galaxy 
catalogue. This number of lens systems was chosen to ob- 
tain a reasonable statistical sample while keeping the time 
needed for computations tractable. The total number of 
'trial' sources was Attempt = 110000 in a total area of 

Aattempt 

= 0.1384deg 2 . Of these 'trial' sources, a number 
unlcns od = 60 ± 7 were unlensed but above the detection 
flux limit Fo. The total number of sources in the total area of 
A survey = 651 deg 2 that would create 400 lenses 4 is 282021, 
yielding a fraction of /i ens = 1.4 ± 0.18 x 10~ 3 of all 'ob- 
served' sources as lensed. This fraction is in good agreement 
with the lens fraction observed in CLASS. 



4-1.1 Lens velocity dispersions and morphologies 

We determine the average velocity dispersion of the 400 
mock lens systems to be a — 164.2±2.6 kms -1 . This value is 
consistent with but slightly lower (by about 5 kms -1 ) than 
what is observed in the CLASS survey. 

The distribution of lens velocity dispersions for 2dF 
lenses expected in a radio selected survey is shown in Fie. 1121 
The distribution for all systems (black) and for 4-image sys- 
tems only (red) is shown. We also show the location of lens 
galaxies in the sample in Fig.|l]in 5 12.3.31 (black dots). Even 
though the lenses lie on the fundamental plane, the complete 
sample of lens systems has a significantly higher velocity 
dispersion than the elliptical galaxy population as a whole. 
However, the velocity dispersion function for 4-image lenses 
follows the distribution of the galaxy velocity dispersion for 
all galaxies much more closely. In the complete sample, a 
total fraction of about 80 ± 5% of lenses are ellipticals with 
a bulge-to-total light ratio of > 0.6. This fraction decreases 
to 70 ± 10% for the 4-image system subsample. The fact that 
these two fractions are very similar suggests that magnifica- 
tion bias is stronger for early-type massive galaxies with a 
moderately asymmetric potential, rather than less massive 



3 The error bar was estimated by running a set of 10 trial runs, 
and noting the value of -/V un i enae( j in each case. 

4 Here, we assume that no sources are lensed that lie outside the 
area Attempt 
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Figure 12. The distribution of lens velocity dispersion in a mock 
lens sample selected in the radio. The black, hatched histogram 
shows the velocity dispersion distribution of all lens systems with 
an average velocity dispersion of a = 164.2 ± 2.6 km s _1 , with the 
majority of lens galaxies (~ 80%) being ellipticals. The red, filled 
histogram shows the statistics for 4 image lenses only. The solid, 
thin black line shows the distribution of velocity dispersions of all 
ellipticals in our mock 2dF catalogue. 
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Figure 13. The expected distribution of image separations in a 
radio selected survey for 2dF lenses. As in the previous figure, the 
black hatched area shows the histogram for all lens systems. The 
red, solid region shows the statistics of 4 image systems only. In 
contrast to Fig. llUl in §[3| this calculation includes magnification 
bias. 

late-type systems with strong ellipticity in the potential due 
to a massive disc. 

4-1.2 Image separations and image geometries 

Using the above selection criteria, we find that the average 
image separation is A# max = 3.3 ± 0.2". Taking into ac- 
count the expected difference in redshift distributions of the 
lenses, z ~ 0.1 and zclass ~ 0.4, and noting that the image 
separations scale as A# max oc Dls/-Dos this compares well 
with the observed image separations in the CLASS sample 

A0CLASS = 1.13". 

The distribution of image separations produced by 2dF 
lenses expected in a radio selected survey is shown in Fia. 1131 



We also show the image separation statistic for a subsample 
consisting of 4-image systems. The average maximum image 
separation decreases to A9 — 2.57T0.6" when only 4-image 
systems are considered. Our model for lenses in the 2dF pre- 
dicts a fraction of 70 ± 5% of double image systems, with a 
fraction of ~ 30 ± 5% being lens systems with four images. 
Depending on the particular choice of known lens systems 
to compare with, these numbers are cons istent /inconsistent 
with o bservations. The sample used by iRusin fc Tegmarkl 
(l200ll) contains seve n quadruples and five doubles, whereas 
the sample used by IChad i2003f) and IChad i2005t) contain 
ten lenses with four quadruple image systems. These sam- 
ples differ in their selection criteria, and since the number 
of known lens systems is low it is difficult to draw any firm 
conclusions. From our analysis, it does appear that the num- 
ber of observed quadruple systems is slightly higher than 
what would be ex pected, even if we compare with the sam- 
ple of lChad ll2005l) . Since our model includes ellipticities and 
disc components it is difficult to explain any further asym- 
metry in the lensing potential needed to increase the frac- 
tion of quadruples with anything other than the presence of 
other mass concentrations in the vicinity of the primary lens. 
These e nviro nmental effects are discussed in more detail in 
iMollerl j200d) . 



5 DISCUSSION AND CONCLUSIONS 

Surveys, like the 2dF and the SDSS are likely to contain a 
large number of gravitational lens systems. We predict the 
number of lenses and lensing statistics expected in the 2dF 
survey. The number of observed lensing systems is a function 
of the statistical properties of the foreground galaxy popu- 
lation, the source population, selection biases and cosmo- 
logical parameters. Taking the cosmological parameters as 
known we calculate the expected lensing statistics for point 
sources at a fixed redshift using realistic multi-component 
galaxy models. Our fiducial galaxy consists of an NFW halo, 
an elliptical bulge and an inclined disc. Detailed galaxy prop- 
erties are predicted using semi-analytic modeling of galaxies 
in the Millennium Run, currently the largest available N- 
body simulation. Using ray-tracing techniques in combina- 
tion with the galaxy model we predict the statistical lensing 
properties of large surveys with realistic galaxy lens models. 

We calculate the lensing cross section for a complete 
sample of 70,000 galaxies at redshifts of zi cns < 0.2 that ful- 
fill the selection criteria of the 2dF survey. With background 
sources at a fixed redshift of z SOUICC — 1 we demonstrate that 
the main contribution to the lensing cross section comes 
from the baryonic component of luminous, early-type galax- 
ies. It is important to note here that this does not mean that 
the lensing cross section is insensitive to the DM component 
in these galaxies. In fact the distribution of galaxy masses 
at a fixed luminosity, n(M v i T \L) for a lens sample is different 
from that for the complete population of early-type galax- 
ies. As shown in Fig.[5]- the mass function of lenses peaks 
at a value that is a factor of about 2 higher than the over- 
all galaxy sample. If magnification bias is not included, we 
find that, for systems with maximal image separations be- 
tween 1 — 3", the disc component of spirals increases the 
lensing cross section by a factor of ~ 10. This shows that 
disc components of spirals are an important contribution to 
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statistical lensing calculations. However, when we do include 
magnification bias we found that this effect is weaker due to 
a trade-off between increase in cross section for late-type 
galaxies on the one hand, and a much stronger magnifica- 
tion bias for four image geometries in moderately elliptical 
massive galaxies on the other (cf. discussion in § 3.3 and 
§4.1.1). 

To make firm predictions with clear selection criteria 
we simulate a sample of 400 radio selected lens systems 
that fulfill the CLASS selection criteria in terms of flux and 
magnification ratios. We find that, in surveys like the 2dF 
or SDSS, including magnification bias, about 1.4 x 10~ 3 of 
all background radio sources are lensed by galaxies within 
the survey. This number is fairly independent of the inclu- 
sion of a disc component, but does depend crucially on the 
presence of the baryonic bulge, and even more strongly, on 
the magnification bias. The number of lenses per observed 
source increases from 1.1 x 10~ 5 without magnification bias 
to 1.4 x 10 -3 including the effect of magnification bias. Us- 
ing our model we also calculate the relative incidence of 4 
and 2 image systems. The inferred fraction of 30 ± 5% of 
quadruple systems is consistent with observed values. 

The key result of our work is that samples of lens galax- 
ies expected in surveys like the 2dF and the SDSS distin- 
guish themselves from the overall galaxy population in sev- 
eral important ways: 

• A fraction of ~ 80% of lenses are expected to be of early- 
type morphologies, compared to only ~ 35% in the overall 
galaxy sample. This number is reduced slightly if only four 
image systems are considered. 

• The average velocity dispersion of lensing galaxies is 
a — 165 km s -1 as opposed to 150kms _1 for all early-type 
galaxies. Furthermore, the distribution is dramatically dif- 
ferent - the distribution for lens galaxies has a maximum 
at a ~ 200kms _1 , with a long tail towards lower velocity 
dispersions. 

• The peak of the luminosity distribution of lens galax- 
ies is a factor of ~ 2 higher than that of the total galaxy 
population. 

• For a given luminosity, lens galaxies reside in haloes 
that are more massive by a factor of ~ 2 on average than 
the overall sample of ellipticals. 

Interestingly, we also find that the properties of the lens 
sample differ depending on whether 4 or 2 image systems 
are considered. Initially surprising is the result that 4-image 
systems tend to have lower velocity dispersions; they trace 
the velocity dispersion distribution of the overall early-type 
galaxy population much better. This can be explained by 
noting that the asymmetries in the potential have a larger 
effect for lower mass galaxies - the ratio of 4: 2 image sys- 
tems is therefore expected to be much higher for galaxies 
with lower velocity dispersions. In addition, the predominant 
cause for 4-image lens geometries is a modest ellipticity in 
the potential, which is more likely to be created by moder- 
ately inclined discs, and hence by later-type galaxies which 
tend to have lower velocity dispersions. Our assumption that 
asymmetries in the potential are due to the baryonic com- 
ponents and that the dark matter haloes are spherical may 
also affect this result. However, since lens galaxies in the cat- 
alogue are dominated strongly by the baryonic component, 
this should not change the statistics significantly. 



All these differences between lens samples and the un- 
derlying galaxy population are very important to under- 
stand if lensing studies are to be used to obtain a better 
understanding of the galaxy evolution. This holds both for 
statistical studies but also for studies of individual systems. 
For example, single cases of lenses for which the mass dis- 
tribution is steeper than the average profile predicted by 
simulations at the same mass do not point to any incon- 
sistency; it is absolutely crucial that the full predicted dis- 
tribution functions are considered. Our statistical analysis 
demonstrates clearly that lensing galaxies are a biased 
population. They preferentially sample the high mass end 
of all galaxies. This systematic, in fact, limits the unbiased 
study of galaxy evolution using a sam ple of lensing galaxies 
jMitchell et alJl2005t lofek et al.ll2003D . 

In our analysis we do not ex amine the effe ct of the en- 
vironment. In a recent study bv lOeuril a2006i) it has been 
shown that the presence of other mass concentrations in the 
vicinity can increase the lensing cross sections for high im- 
age separations by several tens of p ercent. Intere s tingly , this 
result was also previously found bv lTurner et alJ il984f) who 
used a very simple mass-sheet model to calculate the effect 
of environment. We address the effect of environment, and 
in particular the incidence and statistical properties of lens 
galaxies in groups in a subsequent publication. It is worth- 
while to note here, that for low redshifts, the environment 
is expected to have a much smaller effect, as the ratio be- 
tween expected image separations of galaxy lenses and the 
typical distance to nearby galaxies is much smaller at lower 
redshifts (at z ~ 0.05 a distance of 1" corresponds to about 
lkpc, at z ~ 0.5 to about 6kpc). 

Our calculations show that selection effects, and in par- 
ticular magnification biases, have a strong effect on lens 
statistics. Given that magnification bias is important, it does 
appear that the distribution of the sources may also strongly 
influence lensing statistics. For radio selected samples the 
selection criteria can be defined clearly and the source lu- 
minosity function is known to an adequate degree to pre- 
dict lensing statistics reasonably well. In the optical or in- 
frared, however, both selection criteria and source properties 
are not as well understood. To use lenses selected at these 
wavelengths in a statistical way requires careful analysis to 
disentangle the various effects, like source luminosity func- 
tions, magnification bias and source redshift distributions. 
Simulating real imaging data in various wave-bands will be 
needed to test the particular predictions of a given model of 
the lens and source population. 

The basis for this study has been to use one of the best 
models of galaxy formation currently available and test the 
predicted lensing properties of such a sample. Our choice 
of a semi-analytic model was motivated strongly by the ad- 
vantage that the relation between luminous and the DM is 
predicted in a way that is self-consistent, easy to interpret 
and use. It is important to point out, though, that these 
models do have limitations. Even though the predictions at 
low redshift appear to be consistent with observations, it 
is not clear if current semi-analytic models predict all the 
observed properties of galaxies at higher redshifts correctly. 
It is encouraging for semi-analytic models that our results 
are consistent with the current known statistical lens sam- 
ple, but we add the caveat that the known lens sample is 
small. Testing models of galaxy formation and evolution us- 
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ing strong lensing statistics seems a promising approach; 
future lens surveys will increase the number of lens systems 
known by a large number and if the selection criteria in these 
surveys can be understood, much can be learnt about the 
relation between mass and light on galaxy scales. 



6 SUMMARY 

We predict the lensing properties of galaxies in low-redshift 
surveys with z <J 0.2 using realistic galaxy mass models 
based on a recent semi-analytic simulation of galaxies forma- 
tion within the Millennium Run N-body simulation. Using 
ray-tracing techniques we calculate the lensing cross section 
for 2 and 4 image systems for a multi-component lens model 
consisting of a DM halo, a bulge and a disc component. Our 
key results can be summarised as follows: 

• The predicted lensing rate for radio sources in low red- 
shift surveys like the 2dF is 1.4 x 1CT 3 . 

• A fraction of 80 ± 5% of all lens galaxies are predicted 
to be ellipticals. 

• The predicted average maximum image separation is 
3.3 ± 0.2 arcsec. 

• The average velocity dispersion of lenses is a — 164 ± 
2.6kms~\ 

• The baryonic component is the most important contrib- 
utor to the lensing cross sections for galaxy scale masses, 
but the dark matter distribution also affects the lensing 
statistics of composite lens models significantly. 

• Lens galaxies are more luminous and, for a fixed lu- 
minosity, reside in dark matter haloes that have masses a 
factor of 2 higher than the overall galaxy population. 

• The velocity dispersion distribution of lens galaxies is 
shifted to significantly higher values of the velocity disper- 
sion with respect to that of a complete sample of galaxies. 

• Four image systems have, on average, a lower veloc- 
ity dispersion, a later type morphology and a lower average 
maximum image separation than double image systems. 

Our quantitative predictions are entirely consistent with the 
statistics in the CLASS survey. We pioneer the use of de- 
tailed semi-analytic models of galaxies to predict statisti- 
cal lensing properties. This approach provides testable pre- 
dictions for the lensing statistics and conversely, statistical 
lensing can provide useful constraints on galaxy formation 
models in the near future. In ongoing work we address in 
detail the effect of environment for lensing statistics. 
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